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Light quark masses from UKQCD's dynamical simulations with 
0(a) -improved Wilson fermions 

Derek J. Hepburn'', UKQCD Collaboration. 

'^School of Physics, University of Edinburgh, Edinburgh EH9 3JZ, UK. 

I present preliminary results on the light quark masses from a partially quenched analysis of UKQCD's dynam- 
ical datasets. 



1. Introduction 

The determination of the masses of the light 
quarks represents an important challenge in par- 
ticle physics today. They are among the least well 
determined parameters of the Standard Model, 
having uncertainties of the same order as their 
values as given in the particle data book. 

Lattice QCD is in a position to improve this sit- 
uation because it allows one to study the quark 
mass dependence of physical observables such as 
hadron masses. The physical values of these 
observables can then be used to fix the quark 
masses. Unfortunately current simulations can- 
not probe the up/down quark mass range directly 
and so some form of extrapolation is necessary. 
Chiral perturbation theory is at present the best 
procedure for providing a link to the low mass 
regime. The simulations studied in this paper 
however are still too heavy for x^"^ and a simpler 
extrapolation ansatz is used instead. 

In this paper I present results for the average 
up/down quark mass, mi, and the strange quark 
mass, TUs, from a partially quenched analysis of 
pseudoscalar and vector meson masses for each 
of UKQCD's dynamical fermion datasets. The 
suitably averaged values of the hadron masses 
used to fix the quark masses are as follows: 
= 137.3 MeV ; Mk = 495.7 MeV ; M^. = 
893.9 MeV and Mj, = 1020 MeV. 



perturbatively 0(a)-improved Wilson fermion ac- 
tion and standard Wilson glue. Details of these 
simulations have been discussed previously in 
iQJ^. Additionally the set of Kvai values for 
the datasets a.t (3 — 5.2,Ksca — 0.13565 and 
P — 5.2,Ksca — 0.1358 have been increased to 
{0.1342, 0.1350, 0.1355, 0.13565, 0.13580}. Pseu- 
doscalar and vector meson correlation functions 
have been generated with all non-degenerate com- 
binations of these n^ai values. For the Kgoa = 
0.1358 dataset these correlators have also been 
calculated at Kvai — 0.13595. 

3. Hadron masses 

The procedures used here to extract hadron 
masses from correlation functions have been de- 
scribed in detail in ||l|j2) . In brief the masses of the 
pseudoscalar and vector mesons were determined 
from double cosh fits to their corresponding cor- 
relation functions. Local and fuzzed correlators 
were fitted simultaneously to constrain fit param- 
eters and a sliding window analysis used to find 
an appropriate time window for each fit. 

4. Chiral Extrapolations 

The pseudoscalar meson masses for each 
dataset were extrapolated in 1/Kvai using the sim- 
ple ansatz 



2. Simulation Details 

Mps{Ki,K2) = Aps{mi+m2) (1) 
All of the lattice simulations studied here have 111 1 

been carried out on a 16'^ x 32 lattice with a non- ~ ^ps( ^ —) ^ — ) (2) 

Z Ki K2 Kcrit 
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The value of Kcrit was determined as a free pa- 
rameter in this fit. 

For the vector meson the ansatz 

My {mi, 1112) — Av {mi + 7712) + By (3) 

was found to fit the data welL 

An example plot of these chiral extrapolations 
is given in Figure |i| for the Ksca = 0.1358 dataset. 
The lightest pseudoscalar and vector points at 
Kvai = 0.13595 are included in the plot but were 
not included in the extrapolation. The light- 
est vector point clearly deviates from the lin- 
ear behaviour suggesting large finite size effects. 
For more discussion of finite size effects in these 
datasets I refer the reader to [|j . 




Figure 1. Chiral extrapolations of pseudoscalar 
and vector meson masses for the /? — 5.2, Kgca — 
0.1358 dataset. 

5. Fixing the bare light quark masses 

Determination of the light quark masses from 
the results of the chiral extrapolations requires 
setting the scale by some physical quantity. In 
this work the Sommer scale parameter, rg = 
0.49 fm, and Mp = 769.3 MeV were used. 

The average up/down quark mass, mi, was de- 
termined by requiring that the pseudoscalar mass 



equal the physical pion mass at that point i.e. 



Ml,g{mi,mi) = aAU (4) 

The strange quark mass, m^, was then deter- 
mined from each of Mk, AIk' and by solving 

Mpg{mi,ms) = aMx (5) 

Mv{mi,ms) = aMK' (6) 

Mv{ms,ms) = aM^ (7) 



When the lattice spacing was fixed by tq how- 
ever only the Mk method gave a sensible value 
for ms and only these results are given. 

6. Renormalisation and matching 

In this work the values for the renormalised 
quark masses are given in the MS scheme at a 
reference scale of 2 GeV. The relationship be- 
tween the renormalised and bare quark masses is 
given by the relation 

m^{^^) = Z,Y^(/i)(l + b^am,)m, (8) 

where the matching has been performed at a 
scale, fi. Running the masses from the scale, ^, 
to a scale of 2 GeV was achieved using the Math- 
ematica package, RunDec [Q, which implements 
4-loop running from perturbation theory. 

One-loop tadpole-improved perturbative re- 
sults were used for Zm and 6„j reorganised 
in terms of a^'^ (fi) 

Z^^{^,) = [1 ^(8 lnM-15.085)]z.-^(9) 

6,„ = [-1/2 - 0.7850 1 (10) 

The tadpole factor, uq, was taken to be 1/8kc 
here. 

An estimation of the best value for the match- 
ing scale, n = q* 01 , has not been determined for 
this work. Instead the cases = 1/a and fi = n/a 
have been investigated and the variation treated 
as a systematic effect. 

The values of a^{fi) were also determined 
with the RunDec package using values for 
given in 
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7. Results 

Results for mi and nig have been plotted in Fig- 
ures ^ and ^ respectively. Datasets matched in tq 
have been distinguished from the others and a 
matched quenched dataset included for compari- 
son. Only results for /j, = 1/a are plotted as choice 
of matching scale was found to have around a 1% 
effect. The results for ms{M^) are omitted as 
they were almost identical to those for nis {Mk' ) • 

A large dependence on the choice of scale set- 
ting parameter can be seen for both m/ and TOs . It 
should be noted however that a smaller physical 
value for ro would bring the results into better 
agreement. There is also a clear difference be- 
tween ms(M/f) and ms{MK*) though this would 
be expected to vanish in the continuum limit. 

There is some indication of a trend in the 
matched datasets as the sea quark mass is de- 
creased, as expected. When compared to the 
matched quenched point however there is no clear 
signal of unquenching from either quark mass 
though this may improve with statistics. 

The ratio, TOg/m/, was calculated for each 
scheme and found to be between 25.0 and 25.5 
for ms{MK) and around 31.0 for ms{MK')- This 
can be compared with the result from x^'^ 111 of 
24.4(1.5). 
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Figure 2. Average up/down quark masses with 
/I = 1/a. Only statistical errors are shown. 
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Figure 3. Strange quark masses with /i 
Only statistical errors are shown. 
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